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Dissipative phase solitons in semiconductor lasers 
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We experimentally demonstrate the existence of non dispersive solitary waves associated with a 
27r phase rotation in a strongly multimode ring semiconductor laser with coherent forcing. Similarly 
to Bloch domain walls, such structures host a chiral charge. The numerical simulations based on 
a set of effective Maxwell-Bloch equations support the experimental evidence that only one sign 
of chiral charge is stable, which strongly affects the motion of the phase solitons. Furthermore, 
the reduction of the model to a modified Ginzburg Landau equation with forcing demonstrates the 
generality of these phenomena and exposes the impact of the lack of parity symmetry in propagative 
optical systems. 

PACS numbers: Valid PACS appear here 


Dissipative solitary waves have been observed as self- 
localized optical wave packets along the direction of prop¬ 
agation in many out of equilibrium and nonlinear optical 
systems. In spite of their huge variety, many of the obser¬ 
vations reported so far can be cast in two main categories 
depending on the presence or absence of coherent energy 
input, ie the (lack of) phase symmetry of the system 
[T]. In systems with phase symmetry, mode-locked laser 
pulses have been analyzed as dissipative solitons of the 
cubic-quintic Ginzburg Landau equation [2]. Their op¬ 
tical phase can wander in the course of time due to the 
neutral mode created upon the formation of a coherent 
wave. On the contrary, dissipative solitons in forced sys¬ 
tems EH] have been analyzed in the framework of the 
Lugiato-Lefever equation [5], which includes a coherent 
forcing term acting as a phase reference to which solitons 
will lock. In both cases, the use of paradigmatic equa¬ 
tions in addition to system-specific models has allowed to 
formally connect these optical solitary waves to localized 
states as they appear in fluid dynamics, plant ecology, 
granular media or reaction-diffusion systems [6H9]. In 
fact, optical dissipative solitons can often be explained 
as perturbed solitons of the nonlinear Schrodinger equa¬ 
tion (in the weak dissipation limit pTHl2] ) or as locked 
fronts (in strongly dissipative systems [isun]). 

In this Letter, we report on dissipative solitons which 
fundamentally consist of self-confined 2tt phase rotations 
embedded in a homogeneously phase locked background. 
These “phase solitons” are generic features of spatially 
extended oscillatory media under nearly resonant forc¬ 
ing [151 CES] and result from the mismatch between the 
natural periodicity and the forcing. Here, the mismatch 
between the free running laser and the external forc¬ 
ing frequencies leads to the formation of phase kinks 
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FIG. 1: Left: Scheme of the experimental setup. Right: Color 
coded nansi emitted intensity in a space-time representa¬ 
tion. The horizontal axis is a fast time scale normalized to 
the round-trip time of the ring cavity r (close to 3.6 ns) and 
the vertical axis is a slow time scale, in units of r. 


as result of a commensurate-incommensurate transition 
[13 CEB]. This connects our observations with the kink 
solutions observed in many physical systems described 
by the Frenkel-Kontorova model m, such as fluxons in 
Josephson arrays m, local deformations in DNA chain 
[21] , or excitable waves in chemical and biological systems 
[22l 23;. In nonvariational systems, chirality acquired 
through a non equilibrium Ising-Bloch transition implies 
the motion of domain walls [24H26] . A link can thus be 
established with our observations of solitons that indeed 
move in the reference frame of the propagating carrier 
wave (at light speed). Due to the propagative nature 
of the experiment and the non instantaneous semicon¬ 
ductor medium the system lacks parity symmetry, which 
impacts the chiral charge and therefore the motion of the 
phase solitons in this reference frame. 
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The experimental set-up is based on a very large 
(length A ~ 1 m) ring laser featuring a 4 mm-long (/) 
antireflection-coated semiconductor medium inside a free 
space optical cavity (Fig. [I] left) with low mirror trans¬ 
missivity T « 10%. The device operates in a regime 
where the field envelope in the single-mode regime would 
evolve on slower time scale than the active medium [29] . 
In order to avoid transverse effects and directional com¬ 
petition m , the optical cavity includes an aperture and 
an optical isolator. The phase symmetry of the free run¬ 
ning laser is broken by applying an external field provided 
by a grating tunable edge-emitting laser followed by an 
optical amplifier. The emitted intensity is acquired by a 
9.5 GHz photodetector coupled to a 12.5 GHz real-time 
oscilloscope. 

The laser is brought about 10% above threshold and 
emits a superposition of many longitudinal modes. When 
sufficiently strong external forcing is applied, the system 
may lock in phase to the injection beam. In order to 
visualize the dynamics, we construct a co-moving space- 
time representation. We acquire very long time series 
(up to 10 7 points) in real time and then split them into 
segments of equal length, corresponding to the time it 
takes for light to complete a full round-trip of the ring 
cavity. The segments are then stacked on top of each 
other, each of them showing the state of the system at 
a determined round-trip number. In this way, we obtain 
a spatiotemporal representation of the evolution of a 1- 
D system, the horizontal dimension showing the state of 
the system at a given moment and the vertical dimen¬ 
sion its evolution over a discrete time measured in units 
of round-trips. Thanks to very well separated time scales 
(dynamics during one round-trip and evolution over the 
round-trips), the resulting space-time diagrams display 
very clearly the evolution of the system. 

An example is shown on Fig. [l] (right). This regime 
is experimentally attained by choosing initially a very 
strong injected power (sufficient to phase lock the ring 
laser) and setting the detuning to bring the system very 
close to unlocking. Then, upon progressive decrease of 
the injection power and adjustment of the detuning to 
remain close to the unlocking transition, chaotic regions 
spontaneously appear. Here, during the first 1000 round- 
trips, a few spatial regions are in a complex dynamical 
state with only small segments locked to the external 
forcing. Between round-trips 1000 and 3000 a few iso¬ 
lated structures emerge from the chaotic domains and 
drift towards the right, eventually hitting a larger chaotic 
domain which seemingly absorbs them. At rountrip 4000 
a new solitary wave emerges from the chaotic domain 
and drifts first to the left and later to the right, going 
through the boundary on the right and emerging on the 
left about round-trip 7000. The time trace shown on top 
of right panel is a cut through the space-time diagram 
taken at round-trip 6300 (red line). We attribute the 
strong asymmetry between left and right to the lack of 
parity symmetry (also discussed below) in this propaga¬ 
tive system with non-instantaneous medium. 



FIG. 2: Left: Stable phase soliton hosting a 2i r phase slip 
(see inset) at steady state after interaction with a turbulent 
state. Right: Collision and fusion of two solitons. The arrows 
indicate the direction of rotation along the orbit. 


Further information about the nature of these isolated 
structures can be gained by measuring the dynamics of 
the optical phase. This is achieved by measuring simul¬ 
taneously the result of interference between the forcing 
beam and the output beam of the system with three 
different dephasing conditions via a 3 x 3 fiber coupler, 
a simple formula allowing us to compute instantaneous 
intensity and relative phase [31]. This measurement is 
shown in Fig. [2] On the left, the space-time diagram 
shows the coexistence of an isolated structure with a 
chaotic domain, which does not survive the collision with 
the former. The field dynamics is represented in the in¬ 
sets in the (Re(-E), plane. During the first 1000 

round-trips (lower inset), the regions visited by the sys¬ 
tem cover a large part of the plane, indicating that the 
dynamics cannot be well described on this bidimensional 
phase space. On the contrary, at the end of the space- 
time diagram (top inset) most of the system is locked to 
the external forcing and locally circles around the ori¬ 
gin, accumulating 2tt in the field phase. Localized Bloch 
domains carry such a chiral charge but in the present 
case only one stable stationary state exists and the lo¬ 
calized states take the form of a single phase kink [32] . 
As expected [331 , the dynamics of the field in presence of 
solitons projects very neatly on a bidimensional space. 

Contrary to the limit case of a pure phase dynamics 
described by the overdamped Sine-Gordon equation ng, 
the phase rotation is accompanied here by intensity dy¬ 
namics, which takes the trajectory of the system very 
close to the origin, where the chiral charge can disappear 
via the occurrence of a defect m- This provides an in¬ 
terpretation for the right panel of Fig. [2] Two identical 
phase solitons initially propagate at a fixed distance but 
at round-trip 3200, one of them changes speed and ap¬ 
proaches the other until they collide. This interaction 
does not lead to a single kink with a 47r phase rotation 
preserving the total chiral charge. Rather, we observe 
that a single phase soliton with 2tt charge emerges from 
the collision. 

Fig.0 shows the spontaneous nucleation of a phase 
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FIG. 4: Homogeneous synchronous states curve for fi— 1.01, 
a = 3 1 0 = —2.97. The blue trajectory corresponds to the 
phase soliton numerically excited for y — 0.0014 as in Fig. [5] 
The points denoted A, B and C represent a stable node, an 
unstable saddle, and an unstable focus, respectively. 


FIG. 3: Two phase solitons emerge from turbulent regions. 
Each of them hosts a 2tt charge, which is already present in 
the preexisting chaotic domains. 


soliton from a turbulent state. Although the phasor plot 
in the lower inset suggests the presence of multiple tra¬ 
jectories circling around the origin, only the 2tt kink at 
the turbulent domain border survives at the end giving 
rise to a phase soliton. This is actually the typical way 
in which phase solitons appear in our experiment. The 
space-time diagram of Figs. [3] and [l] are strongly reminis¬ 
cent the kink breeding regime documented in [IB] , Yet 
we note that the boundaries of the chaotic domains are 
always markedly different, which we relate to the lack of 
parity symmetry in the spatial dimension. 

To interpret the experimental results, we extend the 
semiconductor models of [34, 35' to include the field lon¬ 
gitudinal propagation in a coherently driven, unidirec¬ 
tional cavity, partially filled by an active medium. The 
phenomenological model for the semiconductor micro¬ 
scopic susceptibility allows for the dependence on fre¬ 
quency and carrier density of the refractive index and 
the gain line. 

Compared to the Haus model for pulse propagation 
in a laser with saturable absorber [36], our approach 
provides a more realistic description of the nonlinear 
radiation-matter interaction in the semiconductor mate¬ 
rial. Thanks to the introduction of material polarization, 
we do not need to introduce an ad hoc spectral filtering 
section to account for material gain dispersion and avoid 
pulses collapse m as happens in lumped-element meth¬ 
ods [38] . 

The model, including the low transmission limit and 
neglecting diffraction, yields a set of equations for the 
normalized polarization P and carrier density D both 
evolving on a faster time scale than the field envelope E, 


which read 

d r E = —rjod^E + a [y — (1 + iO)E + P ], (1) 

d T P = [T(D)(l-ia)p2iS{D)][(l-ia)ED-P}(2) 

d r D = b[fi-D-(E*P + EP*)/ 2], (3) 

with periodic boundary conditions E(0,r) = E(l,r). 
The spatial variable is scaled to the length of the ac¬ 
tive medium, i.e. r] = z/Z, and the temporal variable r 
is actually the retarded time t + z(A — l)/{cl) scaled to 
the polarization decay time r^, the parameter 770 = cr^/A 
is the ratio of the cavity free-spectral range to the gain 
linewidth, which is very small for our long cavity, mean¬ 
ing that even close to threshold a large number of lon¬ 
gitudinal modes may experience gain. The normalized 
decay rates a and b are the polarization-to-photon and 
polarization-to-carrier lifetime ratios, respectively; y is 
the amplitude of the injected field (assumed real), a is 
Henry factor, fi is the scaled pump parameter, fi = 1 
is the free running laser threshold and the value we use 
(fi = 1.01, not critical) has been chosen as the one leading 
to time traces most similar to those of the experiment. 
6 is the detuning between the frequency of the injected 
field (taken as the reference frequency) and the closest 
ring cavity resonance, multiplied by the photon lifetime. 

The functions T(D) = 0.276 + 1.016 D, 5{D) = 
—0.169 + 0.216 D account for the finite gain linewidth 
and for the frequency detuning between gain and cav¬ 
ity, both dependent on the carrier density D and were 
parametrized in previous work by fitting the microscopic 
susceptibility of the active medium [34, [35]. Equa¬ 
tions 0 -@ admit the stationary solution E = \E H 
P = P s , D = D s , with D s = ii/( 1 + X), P s = 
(1 - ia)E s D s , 


4> = arctan [(aD s + 0)/(D s — 1)] , (4) 

y 2 = x[(l-D s ) 2 + (0 + aD s ) 2 ] , (5) 

and X = \E S \ 2 . The input-output relation of Eq. ([5]) 
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FIG. 5: <7 = 3- 1(T 6 , b = 5 • 1(T 4 ,7/ = 0.0014. The other 
parameters are as in Fig. [4] Left: Stable phase soliton seeded 
by a 27r phase slip superimposed to the homogeneous syn¬ 
chronous state A. Bottom left: held intensity profile within 
one roundtrip. Right: Two coexisting phase solitons at steady 
state which display the same phase portrait (insets) 


depends on the parameter /i, <a, 0 and it can be S- 
shaped (three-positive roots, A, B and C in Fig. [4] for 
y = 0.0014), which turns out to be a necessary condition 
for the existence of phase solitons mm- In particular, 
the three fixed points in the subspace (R e(E), Im (E)) 
are, respectively, a stable node (A), a saddle (B) and 
an unstable focus (C). The corresponding phase por¬ 
trait makes the system excitable in absence of propa¬ 
gation ( d/dz = 0) as A and B approach each other 
ns Hang. Linear stability analysis of the input-output 
relation shows that for this parameter set the phase 
locked solution is stable on the whole upper branch, wich 
disappears via a saddle-node bifurcation at the turning 
point, the lower branch and the negative slope branch 
being always unstable. The shaded region in y on Fig. [4] 
corresponds to the stability domain of phase solitons, 
which grows very fast with /i. Below the shaded region 
the locked solution becomes unstable and the system de¬ 
velops an irregular spatiotemporal behavior, while above 
the system uniformly locks to the forcing, preventing the 
existence of phase solitons. 

In the numerical simulations a stable phase soliton can 
be obtained with initial conditions corresponding to the 
locked state A , to which we superimpose a positive phase 
kink of 2ir along z for the field E. The final state is a 
traveling pulse with a 2tt chiral charge and a width of few 
hundreds of ps, which reproduces well the experimental 
results (Fig.[5j left). Its phasor representation in the sub¬ 
space (Re(E), Im(E)) is shown in Fig. [I] by a blue line 
lying in the shaded plane, which confirms that the dy¬ 
namics can be effectively embedded in the 2D subspace. 

Multiple structures separated by varying distances can 
also exist simultaneously (Fig.[5j right). Within paramet¬ 
ric regimes compatible with the experiment, any struc¬ 
ture missing the chiral charge or with the wrong charge 
sign is invariably unstable, though long transients can 
occur. 

It has been shown m that a forced two-level laser 


model (Maxwell-Bloch equations) close to threshold can 
be reduced to a modified Ginzburg-Landau equation 
which admits a free energy equivalent to the potential 
of the Frenkel-Kont or ova model mm- Under certain 
parametric conditions the laser dynamics turns out to be 
slaved to the phase, which obeys an overdamped sine- 
Gordon equation |18|. Although the latter limit does not 
apply to our experiment, Eqs. 0© can be reduced to 
a modified forced Ginzburg-Landau equation using ex¬ 
perimentally meaningful approximations. Precisely, we 
assume that the laser is very close to threshold setting 
fi — 1 = e « 1, and that the injected amplitude y and the 
detuning between the free running laser frequency and 
the injected frequency, which in our model is 9 + cq are 
both of order e. A multiple scale analysis truncated at 
order e 2 has allowed us to eliminate the polarization and 
carrier density while preserving the main features of the 
semiconductor susceptibility (finite width and position of 
the gain peak). The final equation for the field envelope 
E reads 


rjodrjE A- d r E = a {y + [(i - 1 - + 6)] E 

-{l-ia)\E\ 2 E 

+ d [rfcdl - 2iS(l)r] 0 d ri \ E) , (6) 

with d = [r(l) 2 (l — ia)] 1 . The equation differs from 
that of [18] mainly for the last term, which describes a 
nonsymmetric gain line, and for the presence of the Henry 
factor a which implies that the medium is dispersive even 
at resonance, two typical features of semiconductors. In 
a framework moving at velocity 770 (i.e. co-moving with 
the soliton), Eq. © lacks the parity symmetry 77 ——77 
of the standard Ginzburg-Landau equation due to non¬ 
vanishing S( 1). However, unless that parameter is set to 
an unrealistically large value, Eq. © supports both chi¬ 
ral charges contrary to the experiment and the complete 
model. Therefore, we conclude that the main mechanism 
for chiral selection lies in the non-instantaneous medium 
dynamics, which heavily breaks the 77 — 77 symmetry 
in a propagative system. 

In conclusion, we have observed robust solitary opti¬ 
cal structures which propagate in a highly multimode 
semiconductor ring laser with coherent forcing. We 
have measured directly their optical phase and showed 
that they carry a chiral charge whose sign is decisive 
for their stability. Numerical simulations of a physi¬ 
cal model taking into account both the semiconductor 
material susceptibility and the geometry of the exper¬ 
iment allow the interpretation of the experimental ob¬ 
servations in terms of (multistable) phase solitons. In 
addition, multiple scale analysis has allowed to under¬ 
stand these solitons as the elementary excitations of 
oscillatory media under nearly resonant forcing close 
to a commensurate-incommensurate transition. Finally, 
downscaling the whole experiment to monolithic semi¬ 
conductor ring lasers may lead to robust and ultrafast 
phase solitons for phase information encoding and regen- 
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eration in coherent optical communications OHO]. 
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